Whistler Wave Turbulence in Solar Wind Plasma 
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Abstract. Whistler waves are present in solar wind plasma. These waves possess characteristic turbulent fluctuations that are 
characterized typically by the frequency and length scales that are respectively bigger than ion gyro frequency and smaller 
than ion gyro radius. The electron inertial length is an intrinsic length scale in whistler wave turbulence that distinguishably 
divides the high frequency solar wind turbulent spectra into scales smaller and bigger than the electron inertial length. We 
present nonlinear three dimensional, time dependent, fluid simulations of whistler wave turbulence to investigate their role 
in solar wind plasma. Our simulations find that the dispersive whistler modes evolve entirely differently in the two regimes. 
While the dispersive whistler wave effects are stronger in the large scale regime, they do not influence the spectral cascades 
which are describable by a Kolmogorov-like fc~ 7 / 3 spectrum. By contrast, the small scale turbulent fluctuations exhibit a 
Navier-Stokes like evolution where characteristic turbulent eddies exhibit a typical fc -5 / 3 hydrodynamic turbulent spectrum. 
By virtue of equipartition between the wave velocity and magnetic fields, we quantify the role of whistler waves in the solar 
wind plasma fluctuations. 
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1. INTRODUCTION 

The solar wind is an excellent in-situ laboratory for in- 
vestigating nonlinear and turbulent processes in a mag- 
netized plasma fluid since it comprises a multitude of 
spatial and temporal length-scales associated with an ad- 
mixture of waves, fluctuations, structures and nonlinear 
turbulent interactions. The in-situ spacecraft measure- 
ments 01 0] reveal that the solar wind fluctuations, ex- 
tending over several orders of magnitude in frequency 
and wavenumber, describe the power spectral density 
(PSD) spectrum that can be divided into three distinct 
regions J2, Si- The frequencies, for instance, smaller 
than 10 5 Hz lead to a PSD that has a spectral slope 
of -1 . This follows the region that extends from 10 5 
Hz to or less than ion/proton gyrofrequency where the 
spectral slope exhibits an index of -3/2 or -5/3. Smaller 
than ion gyro radius (fep, ^> 1) and temporal scales big- 
ger than ion cyclotron frequency co > (0 C i = eBo/m e c, 
(where k,pi,(0 C i,e,Bo,m e ,c are respectively characteris- 
tic mode, ion gyroradius, ion cyclotron frequency, elec- 
tronic charge, mean magnetic field, mass of electron and 
speed of light) exhibit a spectral break where the iner- 
tial range slope of the solar wind turbulent fluctuations 
varies between -2 and -5 J2, [H. The onset of the sec- 
ond or the kinetic Alfven inertial range is still elusive to 
our understanding of the solar wind turbulence and many 
other nonlinear interactions. Specifically, the mechanism 



leading to the spectral break has been thought to be ei- 
ther mediated by the kinetic Alfven waves (KAWs) @], 
or by electromagnetic whistler fluctuations JUS], or by 
a class of fluctuations that can be dealt within the frame- 
work of the HMHD plasma model U [k^]. Stawicki 
et al |7[] argue that Alfven fluctuations are suppressed 
by proton cyclotron damping at intermediate wavenum- 
bers so the observed power spectra are likely to consist 
of weakly damped magnetosonic and/or whistler waves 
which are dispersive unlike Alfven waves. Moreover, tur- 
bulent fluctuations corresponding to the high frequency 
and kpj ^> 1 regime lead to a decoupling of electron mo- 
tion from that of ion such that the latter becomes unmag- 
netized and can be treated as an immobile neutralizing 
background fluid. While whistler waves typically survive 
in the higher frequency (and the corresponding smaller 
length scales) part of the solar wind plasma spectrum, 
their role in influencing the inertial range turbulent spec- 
tral cascades is still debated lO[H[3[Il[I!. 

In this paper, we focus on understanding the nonlinear 
turbulent cascades mediated by whistler waves in a fully 
three dimensional geometry. Our objective is to inves- 
tigate the role of whistlers in establishing the turbulent 
equipartition amongst the modes that are responsible for 
the nonlinear mode coupling interactions which critically 
determine the inertial range power spectra. 
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2. WHISTLER WAVE MODEL 

Whistler modes are excited in the solar wind plasma 
when the characteristic plasma fluctuations propagate 
along a mean or background magnetic field with fre- 
quency co > CO C i and the length scales are c/(o p i < £ < 
c/cOp e , where (Opi,CO pe are the plasma ion and elec- 
tron frequencies. The electron dynamics plays a criti- 
cal role in determining the nonlinear interactions while 
the ions merely provide a stationary neutralizing back- 
ground against fast moving electrons and behave as scat- 
tering centers. The whistler wave turbulence can be de- 
scribed by the electron magnetohydrodynamics (EMHD) 
model of plasma 11611 . The three-dimensional equation 
of EMHD describing the evolution of the magnetic field 
fluctuations in whistler wave, 

4-(B-^ 2 V 2 B)+V,-V(B-^ 2 V 2 B)- (1) 
at 

(B - c/ 2 V 2 B) • VV e = lld*V 2 B. 

The length scales in Eq. (Q]) are normalized by the elec- 
tron skin depth d e = cj co pe i.e. the electron inertial length 
scale, the magnetic field by a typical amplitude Bo, and 
time by the corresponding electron gyro-frequency. In 
Eq. (Q}, the diffusion operator on the right hand side 
is raised to 2n. Here n is an integer and can take n — 
1,2,3, The case n — 1 stands for normal diffusion, 
while n = 2,3, • • • corresponds to hyper- and other higher 
order diffusion terms. 

The linearization of Eq. (Q~|) about a constant magnetic 
field Bq yields the dispersion relation for the whistlers, 
the normal mode of oscillation in the EMHD frequency 
regime, and is given by 

where co co = eBo/mc and k 2 = k 2 + k 2 . From Eq. ([TJ, 
it appears that there exists an intrinsic length scale d e , 
the electron inertial skin depth, which divides the entire 
spectrum into two regions; namely short scale (kd e > 1) 
and long scale (kd e < 1) regimes. In the regime kd e < 1, 
the linear frequency of whistlers is co^ ~ kyk and the 
waves are dispersive. Conversely, dispersion is weak in 
the other regime kd e > 1 since C0k ~ k y /k and hence the 
whistler wave packets interact more like the eddies of 
hydrodynamical fluids. 

3. SIMULATIONS 



whistler wave turbulence code that we have developed 
at Center for Space Plasma and Aeronomic Research 
(CSPAR), the University of Alabama in Huntsville 
(UAH). Our code numerically integrates Eq. (Q~|). The 
spatial descritization employs a pseudospectral algo- 
rithm 
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based on a Fourier harmonic expan- 
sion of the bases for physical variables (i.e. the magnetic 
field, velocity), whereas the temporal integration uses a 
Runge Kutta (RK) 4th order method. The boundary con- 
ditions are periodic along the x,y and z directions in the 
local rectangular region of the solar wind plasma. 

Electron whistler fluid fluctuations, in the presence of 
a constant background magnetic field, evolve by virtue 
of nonlinear interactions in which larger eddies transfer 
their energy to smaller ones through a forward cascade. 
According to 11911 . the cascades of spectral energy occur 
purely amongst the neighboring Fourier modes (i.e. lo- 
cal interaction) until the energy in the smallest turbulent 
eddies is finally dissipated gradually due to the finite dis- 
sipation. This leads to a damping of small scale motions. 
By contrast, the large-scales and the inertial range turbu- 
lent fluctuations remain unaffected by direct dissipation 
of the smaller scales. Since there is no mechanism that 
drives turbulence at the larger scales in our model, the 
large-scale energy simply migrates towards the smaller 
scales by virtue of nonlinear cascades in the inertial 
range and is dissipated at the smallest turbulent length- 
scales. The spectral transfer of turbulent energy in the 
neighboring Fourier modes in whistler wave turbulence 
follows a Kolmogorov phenomenology lfl9l EH H3l that 
leads to Kolmogorov-like energy spectra. We find from 
our 3D simulations that whistler wave turbulence in the 
kd e < 1 and kd e > 1 regimes exhibits respectively fc~ 7 73 
and £~ 5 / 3 (see Fig 1) spectra. The inertial range turbu- 
lent spectra obtained from our 3D simulations are fur- 
ther consistent with 2D work iflll l2~ll I22I1 . Interestingly, 
the wave effects dominate in the large scale, i.e. kd e < 1, 
regime where the inertial range turbulent spectrum de- 
pletes a Kolmogorov-like &~ 7 / 3 spectrum. On the other 
hand, turbulent fluctuations in the smaller scale (kd e > 1 ) 
regime behave like non magnetic eddies of hydrody- 
namic fluid and yield a A: -5 / 3 spectrum. The wave effect 
is weak, or negligibly small, in the latter. Hence the non- 
linear cascades are determined essentially by the hydro- 
dynamic like interactions. The observed whistler wave 
turbulence spectra in the kd e < 1 and kd e > 1 regimes 
(Figs 1) can be followed from the Kolmogorov-like argu- 
ments 1 19, that describe the inertial range spectral 
cascades. We elaborate on these arguments to explain our 
simulation results of Fig. (1) in the following section. 



Turbulent interactions mediated by the coupling of 
whistler waves and inertial range fluctuations are stud- 
ied in three dimensions (3D) based on a nonlinear 3D 
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FIGURE 1. (LEFT) 3D simulation of whistler wave turbulence in the kd e < 1 regime exhibits a Kolmogorov-like inertial range 
power spectrum close to kT 1 ^ . The simulation parameters are: Box size is L x x L y x L, = 2k x 2k x 2k, numerical resolution is 
N x x N y xN : = 200 x 200 x 200, electron skin depth is d e = 0.015, magnitude of constant magnetic field is fi = 0.5. (RIGHT) the 
kd e > 1 regime depicts a Kolmogorov-like k~ 5 ^ spectrum. d e = 0.15. 



4. WHISTLER WAVE SPECTRA 

The exact spectral indices corresponding to the whistler 
wave turbulent spectra, described by the ideal elec- 
tron magnetohydrodynamic invariant, can be understood 
from the Kolmogorov's dimensional arguments [19, 

Hi 

In the underlying whistler wave model of magnetized 
plasma turbulence, the inertial range eddy velocity is 
characterized typically by v e ~ V x B. Thus the typical 
velocity of the magnetic field eddy B( with a scale size 
t can be represented by v e ~ B( / 1, The eddy turn-over 
time is then given by 

£ £ 2 

T ~ — ~ — . 

V e B e 

This is the time scale that predominantly leads to the 
nonlinear spectral transfer of energy in fully developed 
whistler wave turbulence. 

In the regime where characteristic length scales are 
bigger than the electron skin depth (kd e < 1), the inertial 
range whistler turbulent energy is dominated by the large 
scale fluctuations. The total energy corresponding to the 
turbulent fluctuations in this regime is then given as, 

£~ |B| 2 ~Bf ~v 2 £ 2 . 

The B( represent magnetic field associated with the mag- 
netic field eddy of length I. The second similarity fol- 
lows from the assumption of an equipartition of energy 
in the magnetic and velocity field components of whistler 
waves. The process of equipartition origintes from the 
correlation between the velocity and magnetic field fluc- 
tuations v t ~ k x B, where k = k x x + k y y + k z z is a three 
dimensional wave vector. The latter is further consistent 



with the electron flow speed, in combination with the 
wave perturbed magnetic field, that is used to derive the 
dynamical equation of whistler wave turbulence, i.e. Eq. 
(Q]). This velocity-magnetic field correlation essentially 
produces the velocity field fluctuations that are normal 
to the magnetic field in a whistler wave packet. Conse- 
quently, the energy associated with the velocity and mag- 
netic field for each characteristic turbulent mode evolves 
toward a relationship that satisfies v 2 ~ k 2 B 2 . To quan- 
tify our arguments, we follow the evolution of turbulent 
equipartion in our simulations by computing the follow- 
ing quantity, 

^(O^LdVe^OI 2 -* 2 ^.*)! 2 ). ( 2 ) 
k 

Interestingly, our 3D simulations, describing the equipar- 
tition between the velocity and magnetic field fluctua- 
tions, are consistent with the 2D counterpart. It thus ap- 
pears that the turbulent equipartition is a robust feature of 
whistler waves that is preserved in both 2D and 3D non- 
linear mode coupling interactions. The spectral cascades 
of inertial range turbulent energy is nonetheless deter- 
mined by the energy cascade per unit nonlinear time as 
follows, 

E B\ 
p ~ — ~ — - 

T " I 2 ' 

On assuming that the spectral energy cascade is local in 
the wavenumber space SIM Si, the energy spectrum 
per unit mode yields ~ E a kP. On substituting the en- 
ergy and energy dissipation rates and equating the pow- 
ers of B( and £, we obtain a = 2/3 and j3 = —7/3. This, 
in the kd e < 1 regime, leads to the following expression 
for the energy spectrum E% ~ £ 2 / 3 £~ 7 / 3 . 
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FIGURE 2. Turbulent equipartion between the velocity and 
magnetic fields is observed in our 3D simulations. When the 
characteristic turbulent modes evolve towards equipartition, the 
relationship |v e (k,f)| 2 ~ k 2 \B(k,t)\ 2 is obeyed. Consequently, 
F ■ — > 10~ 7 

5. CONCLUSIONS 

Three dimensional simulations of turbulent cascades in 
solar wind plasma are carried out to quantify the role 
of whistler waves corresponding to the inertial range 
fluctuations that possess characteristic frequency bigger 
than the ion gyro frequency (co > C0 C i) and length scales 
smaller than the ion gyro radius (kpi > 1). In this regime, 
the solar wind plasma fluctuations comprise of unmag- 
netized ions, hence the entire dynamics is governed by 
the electron fluid motions. The rotational magnetic field 
fluctuations in the presence of a background magnetic 
field lead to propagation of dispersive whistler waves in 
which the wave magnetic and velocity fields are strongly 
correlated through the equipartition (v 2 ~ k 2 B 2 ). The lat- 
ter is employed in our simulations to quantify the role 
of whistler waves that are ubiquitously present in the in- 
ertial range in the high frequency (co > ©„•) solar wind 
plasma. Interestingly we find that despite strong wave ac- 
tivity in the inertial range, whistler waves do not influ- 
ence the inertial range turbulent spectra. Consequently, 
the turbulent fluctuations in the inertial range are de- 
scribed by Kolmogorov-like phenomenology IU9I1 . 

It is to be noted that as long as the cascade of energy 
is concerned, kinetic JH 0] and fluid 
processes lead to a similar power law (i.e. E\ ~ fc~ 5 / 3 ) 
in the kd e > 1 regime. This is because the energy cas- 
cade is determined entirely by the convective time scales 
that are associated with the nonlinear term V e • VV 2 B in 
the electron fluid momentum equation. The breakdown 
of fluid-like behavior occurs for the characteristic scales 



kd e 3> l,kp e ~ 1, where p e is electron gyro radius. The 
major difference in the fluid and kinetic models, how- 
ever, arises from wave-particle resonances (or wave par- 
ticle interactions) which are a fully kinetic effect and it is 
beyond the capability of the fluid theory. Since the energy 
spectra are not critically dependent on the wave-particle 
resonances, our fluid model yields spectral laws similar 
to the kinetic model. 

To conclude, consistent with the Kolmogorov-like di- 
mensional argument [19], we find that turbulent spectra 
in the kd e < 1 and kd e > 1 regimes are described respec- 
tively by £~ 7//3 and £~ 5//3 . Our results are important par- 
ticularly in understanding turbulent cascade correspond- 
ing to the high frequency (co > CO C i) solar wind plasma 
where characteristic fluctuations are comparable to the 
electron inertial skin depth. 
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